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Abstract 

We investigate timelike junctions (with surface layer) between spherically symmetric solutions of the 
Einstein-field equation. In contrast to previous investigations this is done in a coordinate system in 
which the junction surface motion is absorbed in the metric, while all coordinates are continuous at the 
junction surface. 

The evolution equations for all relevant quantities are derived. We discuss the no-surface layer case 
(boundary surface) and study the behaviour for small surface energies. It is shown that one should expect 
cases in which the speed of light is reached within a finite proper time. 

We carefully discuss necessary and sufficient conditions for a possible matching of spherically sym- 
metric sections. 

For timelike junctions between spherically symmetric space-time sections we show explicitly that the 
time component of the Lanczos equation always reduces to an identity (independently of the surface 
equation of state). 

The results are applied to the matching of FLRW models. We discuss 'vacuum bubbles' and closed- 
open junctions in detail. As illustrations several numerical integration results are presented, some of 
them indicate that the junction surface can reach the speed of light within a finite time. 

1 Introduction 

Recent measurements of the microwave background radiation support the idea that our universe is highly 
isotropic and homogeneous Cosmological models with these properties are uniquely represented by the 
class of Friedmann-Lemaitre-Robertson- Walker (FLRW) models. The observed flatness of the universe can 
then be explained by employing an inflationary model which suggests an exponential expansion of the 
early universe driven by a scalar field - the inflaton field. 

Nevertheless, it is often speculated that this might only be the local geometry, while over a larger scale 
the universe is inhomogeneous and anisotropic, i.e., the matter content and geometry vary. In particular, it 
appears as if the parameters necessary for life are highly fine-tuned and in order to solve this "fine-tuning 
problem" it was suggested that we live in one of many different FLRW regions - most of them might be 
unsuitable for life. The most prominent example is Linde's Chaotic Inflation Scenario OE], in which the 
different FLRW regions originate from different almost homogeneous Planck-sized regions which experience 
a period of exponential expansion. 

When such models are discussed it is usually assumed that the transition region between two almost 
FLRW regions is very small and can be approximated by a timelike junction hypersurface, the so-called 
thin bubble wall. To find the motion of this hypersurface one has to find the matching surface to the two 
solutions of the Einstein-field equation representing the space-time on each side. 

The matching conditions are of two different types. On the one hand there is the purely geometric 
necessity that 'things fit together' — distances on the junction surface should have the same length when 
measured 'inside' or 'outside'. On a mathematical level this reduces to a matching of the tangential metric 
components. 



On the other hand there are the matching conditions which result from the assumed vahdity of certain 
physical laws, in particular the energy-momentum conservation across the junction surface and the validity 
of the Einstein-field equation on each side. These conditions have been evaluated by C. Lanczos 0, R- 
Dautcourt [Sj , and in a ground-breaking work by W. Israel Q . 

While these equations are in principle valid for any matching of two space-times satisfying the Einstein- 
field equation, it is in practice impossible to handle their complexity except for highly symmetric cases and 
for all practical applications spherical symmetry is assumed. There is a vast amount of literature, see e.g., 

[ziiHiiaiini. 

A good introduction to the standard method for dealing with spherically symmetric junctions is given 
by K. Lake in 10 , where also different matchings, including cosmological voids and vacuum bubbles, are 
studied numerically. In j9j V. A. Berezin, V. A. Kuzmin and I. I. Tkachev considered the generic spherically 
symmetric case and expressed the relevant quantities with invariants. Junctions arising from phase transitions 
in the early universe, including junctions between sections of FLRW sections are discussed. 

Our aim here is to present a new approach to junctions between spherically symmetric space-times and 
to apply the formalism to junctions between FLRW models — analytically and numerically. 

The approach will focus on the geometrical quantities describing the situation, i.e., the distance of the 
junction surface from the centre of symmetry, and not alone on the junction surface radius. In contrast 
to most other studies, we do not evaluate the junction conditions (in particular the Lanczos equation) in 
the original coordinate system or in Gaussian normal coordinates based on the junction surface. Instead 
we introduce new coordinates such that the junction surface is at a fixed (new) "radial" coordinate and all 
coordinates are continuous at the junction surface. The motion of the junction surface is now absorbed into 
the metric components. 

In spherically symmetric cases the Lanczos equation has two non-trivial independent components - an 
angular and a time component. While the first one leads to the well-known evolution equations, there 
seems to be uncertainty about the interpretation of the time component, which is a second order (in time) 
differential equation for the junction surface motion. It has been known that for certain particular cases 
this equation reduces to an identity JI]. Nevertheless, other authors^ suggested that this equation acts as a 
surface equation of state [^j, i.e., determines the surface pressure. Using the presented approach we will show 
that the time component of the Lanczos equation is in fact an identity for all junctions between spherically 
symmetric solutions of the Einstein-field equation. 

It should be pointed out that there are special cases of junctions which could be examined without 
employing junction conditions. If the 7-equation of state and the cosmological constant have no discontinuity 
at the junction surface then the spherically symmetric space-time can be described in terms of the Lemaitre- 
Tolman model. However, the really intersting question is how the junction behaves if the inside and outside 
region have different dynamical behaviour, i.e., different equation of state and cosmological constant. For 
these cases one cannot avoid the use of junction conditions and all numerical examples given in this paper 
will be of this kind. 

Our approach differs in the use of a coordinate transformation to minimize the number of discontinuous 
quantities and to evaluate the standard junction conditions in a more convinient form (leading to H23|l and 
(|24ll ). The main results do agree with well known results in the literature — for example the evolution 
equation for the angular metric component (|30|1 can be found in a similar (though I believe less convinient) 
form in ^U] and the results for vacuum bubbles agree with findings in ^ and |12j . 

This paper will be structured as follows: In section|21we re-examine junction conditions for the matching 
of generic spherically symmetric sections. We will re-derive the evolution and constraint equations using a 
new approach, based on a coordinate transformation which makes all coordinates continuous at the junction 
surface. We investigate the behaviour for small values of the surface-energy density and discuss the special 
case of vanishing surface-energy density. In section|3|we pay particular attention to constraints on the surface- 
energy density and in section^lwe show that the time-component of the Lanczos equation is an identity. This 
is followed by an application to the matching of FLRW sections in section |S1 and three numerical examples 
in section 

1.1 Notation 

We use standard notation and metric signature ( — I — h+). Coordinates are labeled with greek indices, running 
from to 3, where represents the time coordinate. Latin indices label coordinates on the three-dimensional 

^In particular this was suggested for the junction between FLRW models. 



2 



junction hypersurface. As will be seen below, with our choice of coordinates the junction surface is located 
at a fixed radial coordinate R = 1, and hence latin indices take the values 0,2, and 3 (or equivalently t,9, 
and (/))■ 

When considering quantities defined for certain hypersurfaces (like the junction surface) it will be con- 
venient to have a covariant derivative for this subspace. We will use a vertical bar (like in Kab\c) to refer 
to this covariant hypersurface derivative, which is evaluated in the same way as above covariant derivative, 
but with all quantities replaced by the corresponding hypersurface quantities. Any quantity referring to the 
three-dimensional hypersurface of the junction surface will have a superscript '(3)', e.g., ^^^R and ^^^gab- 

We will use a subscript — or -|- sign to indicate whether quantities refer to the in- or outside region, 
respectively. 



2 Matching of generic 0(3)-symmetric sections 
2.1 The coordinate system 

Any spherically symmetric space-time (i.e., having 0(3) symmetry) allows coordinates such that the metric 
takes the form 

ds2 ^ -m^ir, r)dT^ + P{t, r){dr^ + f{T, r)dn^}, (1) 

where *H(t, r) is the so-called lapse function, and dil^ — d9^ + sin^((?)d0^ the line-element on the two- 
dimensional unit-sphere. 

We are interested in the matching of two spherically symmetric space-times, each having a metric of the 
form Generally, the coordinates will not match up at the junction surface and the manifold is described 
by two different coordinate charts - one for the inside region (subscript — ) and one for the outside region 
(subscript +). At any coordinate time (in the inside or outside region) the junction surface itself is assumed 
to be a two-sphere which is described by its coordinate radius in the inside and outside region, a_(T_) and 
a+(T+). 

Since the inside and outside regions are originally described by a metric in the form |^ we will use the 
convention that a dot /prime refers to the proper time/radial derivative with respect to the metric at the 
junction surface, e.g.. 



,• dcf 1 9 , , - 



and l',=^-^l+{T,r) 
L+ or 



r=Q+(T+) 



In order to describe the motion of the junction surface one usually tries to find the evolution of the 
junction surface radius in each coordinate system. Here we want to suggest a different approach: we 
introduce a new coordinate system, such that all the coordinates are continuous at the junction surface 
while only the transverse metric components are discontinuous at the junction surface. The junction surface 
motion is now described by the evolution of the metric components. 

This, however, should not be confused with Gaussian normal coordinates (which are widely used in the 
discussion of junctions) as our coordinate curves =const. will be in general non-orhogonal at the junction. 



Constructing a continuous time coordinate Let us assume that the inside and outside spaces are 
given in terms of their metrics, which take the form Generally the time coordinates for the inside and 
outside region will not match up at the junction hypersurface. 

Since the junction surface is timelike, each value of the time-coordiante (in the inside or outside re- 
gion) identifies a unique spherically symmetric hypersurface of the junction. This establishes a strictly 
monotonically increasing one-to-one relation between the times on each side t+ = F{t-). Setting dr_ == 
dT+ = F'{t)dt introduces the new continuous 'global time coordinate' t. This rescaling is independent of r 
and leaves the form of the metric Q invariant. However, if the original lapse function was constant then 
this re-scaling results in a new time dependent lapse function which contains information about the junction 
surface motion. For example, this will be the case for the matching of FLRW models. 

If on the other hand the original lapse function 01(t, r) , where r is the original time coordinate, depends on 
the radial coordinate then we can write the new lapse function (which makes the time coordinate continuous) 
as 

N{t,r)=F'{m{T,T), 
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Old Coordinates 



New Coordinates 




Figure 1: An illustration of the coordinate transformation. The original radial coordinates r_|_ and r_ are 
rescaled such that in the new coordinate system the junction surface is at a fixed 'radial' coordinate R = 1. 
The time-coordinates are rescaled such that they match up at the junction surface. 



where r — J^F'{t')dt' and t is the new time coordinate. It is now ri{t) which contains information about 
the junction surface motion, while ^ contains information about the background space. In particular, the 
quantity 

N'{t,r) _ W(r, r) 
N{t,r) ~ m{T, r) 

(taken as a function of proper time) does not depend on F' . We will keep these issues in mind when we use 
the lapse function in the following calculations. 

Finally it should be pointed out that there remains the 'gauge freedom' to rescale the new global time- 
coordinate which will be used later. 



Constructing a continuous radial coordinate Now we want to construct new radial coordinates such 
that the junction surface is at a fixed radial coordinate i? = 1. This can be achieved by setting r± = a±(t)R, 
where the subscript -I- refers to outside (i? > 1) and — to the inside region, and a±{t) is the coordinate 
radius of the junction surface at time t. The relation between the old and the new coordinates is illustrated 
in figure H With this new radial coordinate the metrics for the inside and outside region take the form 

ds^ = ~N^{1 - R'^a^fjdt^ + 2aNafRdR dt + l^{a^dR^ + f{aR)dn'^}, (2) 
where a dot indicates the proper time derivative along paths of constant r, 9, 0, i.e., a = jr^tCe. 



2.2 Geometric matching conditions 

At the matching surface R—1 the tangential metric components must be continuous. This gives us the two 
matching conditions 

[If] = 0. (3) 

and 

[N^ia'^l^ ^l)]=0, (4) 

where 

[giR)] lim - lim {g{R)) = lim (g+(r+)) - lim (.g-(r_)). 

These relations identify two quantities which are continuous across the junction surface and we define 

kit,R) N^l-a^PR^ , k = k(t,l)=N^l-aH^ (5) 

and 

w'^lf, (6) 
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which are the tangential metric components. Note that k{t, R) becomes complex for large cPPB? . However, 
this is not of relevance for the problem at hand, since we are only interested in the behaviour around the 
junction surface at i? = 1 where R) is real for timelike junction surfaces. 

It follows from Q that if the surface appears from one side as timelike, it will so from the other side. 
From now on let us assume that the junction surface is a timelike surface, i.e., a±l± < 1. 

There is a remaining gauge freedom: we can rescale the global time coordinate, i.e., multiply the lapse 
functions with a time dependent factor. One particularly useful choice is to rescale the time such that the 
tangential metric component in the timelike direction parallel to the junction surface becomes unity, i.e.. 



k = Ny/l-a^f 



= 1- (7) 

R=l 



To maintain generality we will not assume this choice until explicitly stated (in section r2.4|l . 

The two conditions ^ and Q are of pure geometric character - they have to be satisfied independently 
of the evolution equations at all times. Taking the total derivative of ^ with respect to coordinate time we 
obtain the corresponding restriction on the junction surface motion 



^ = N{{lfr + {IfYal} 



dw 
'dt 



0. (8) 



It should be noted that here N is not independent, but depends via on the junction surface motion. 
Equations and JH)) generally have two solutions for the junction surface motion in terms of the surface 
radius evolution. 

Let us note here that the metric of the timelike hypersurface representing the junction surface is given 

by 

'^^''g^^dx^'dx" = ds| = -k^df + w^dn"^. 

The intrinsic geometry of the junction hypersurface is completely defined by k and w. Nevertheless, w might 
not uniquely identify the position of the junction surface. Only if l+f+ and l-f- are invertible functions of 
a± then the position of the junction surface is indirectly given by the value of w. 

2.3 Lanczos equation and Israel junction conditions 

Let us split the energy- momentum tensor in a regular and a 5- function part, so that 

Tf_,„ 5{iii)Si_,^ + ff^^, (9) 

where T^i, contains the regular part and 77 is a function of the coordinates which vanishes on the junction 
surface, is non-zero everywhere else, and on the junction surface its gradient is a unit vector. The tensor 
S^i, is called the surface stress-energy (or energy-momentum) tensor. The i5-function restricts its influence 
to the junction surface and we assume that it only depends on coordinates on the junction surface, i.e., in 
our case this tensor does not depend on R. The Lanczos equation [Sj relates the surface energy-momentum 
tensor Sf^i, to the jump in the extrinsic curvature K^^, of the junction surface by 

kS'"' = ^^'>g'"'[K]~[K'"'], (10) 

or equivalently (after taking the trace and substituting back for K) 

-[K'^n = -S''" + -^^'>g'"'S, (11) 
K 2 

where K = K'^fj,, S 5''^^, and k = SttG. These two equations imply that the presence of a surface layer 
is equivalent to a jump in the extrinsic curvature, i.e., 7^'' = [X^"] 7^ 0. 

To relate this conditions to the metric components on both sides of the junction surface we have to find 
expressions for the extrinsic curvature. We start with the normal to the junction surface, which is given by 



^ VI - R^a^P 



(12) 

R=l 
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and the unique timelike unit-vector tangential to the junction surface and orthogonal to the spherical sym- 
metric subspace, which is given by 

= i-^f. (13) 

The extrinsic curavure i^^^ = ni^x-j^^h^ ^h'^ where h^^'^ — g^'^ — n^^n^ is the projection tensor onto the 
junction surface, has the two independent components 

„ „ 1 N + 2N'la + aH^N(i/l) , , 

K%^^{al{lfy + {lfy), (15) 
kw 

where we used the coordinate time derivative of to simplify the first expression. 

This form of the extrinsic curvature implies that the surface energy-momentum tensor 5'^'' is diagonal 
and of perfect- fluid form (in the three-dimensional hypersurface space). We introduce the surface energy 
density ps and pressure such that 

S'^^ = (p, +p,)u'^u'^ (16) 
The Lanczos equations l|ll|l are now given by the two equations 

[N{al{fiy + = -^Pswk (17) 

N + 2N'la + d^PN{i/l) = K Qp, + k^. (18) 

The Gauss-Codazzi equations relate the curvature at one point to the extrinsic and intrinsic curvature 
of a hypersurface which passes through this point by 



R^i.XKea^eb^ec^ed" = Rabcd + e{KadKbc - KacK, 



bd) 

and Rf_tiy\K'n^ eb'^ ed'^ — Kbc\d ~ Kbd\c: where is a coordinate basis of the junction surface, a vertical 
bar denotes the covariant derivative with respect to the induced hypersurface metric — g^^, — n^rii^, 
and e = n^n^ equals -1-1 for timelike hypersurfaces, and —1 for spacelike hypersurfaces. For a timelike 
hypersurface these equations lead to (after substituting ((TT|l ') 

[G^^ea^n''] = -~KSa\b (19) 

[G^,.n^'n''] = \[K^ - K^^K^^'^] - kS^^^K^,, (20) 

where K^^ i(limfl._^i+ K^i, + limjj.^i- K^y) and G^jy = — ^9fj.uR is the Einstein tensor. 

The second of these equations represents nothing else than the definition of the surface stress-energy 
tensor S'^'^ (which was substituted) and hence is redundant. Using H16(l to evaluate H19(l we find that the 
only non-vanishing component is the time component, which is given by (using cq^ = it^) 

ku'^p,.^ + kip, + p,)u^\a = ti-'lG^yU^n-'] = [T^.u^n''], (21) 

where the last equality follows from the Einstein equation and = u^(2 ln(Z/));^. Together with an 
equation of state for the surface energy and pressure densities this equation describes the evolution of the 
'matter' on the surface. 



2.4 Matching with surface-layer 



The second geometric matching condition shows that k = ~ a'^P is continuous across the junction 

surface. Hence we can express the junction surface motion in terms of the lapse function by 




«^ = ±\1-U7 ■ (22) 
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For convenience we will choose now k = 1, i.e., the coordinate time corresponds to proper time along 
the curves R = 1,0, </> constant on the junction surface. Setting u N/k and j± sgn(d±) the angular 
component of the extrinsic curvature (|15|l on both sides of the junction surface and equation ||HJ), the derivative 
of the jimction surface radius with respect to coordinate time, are given by 



K% = (W + jwy/u^ - (23) 



dw 



L = — = UW + jw' \/v? -\ (24) 
At 

A valid matching between two spherically symmetric sections, each satisfying the Einstein-field equations, 
must satisfy the two geometric matching conditions (lj3l and I^J) and the two independent components of the 
Lanczos equation f ljlTfl and (|18|l '). We first note, that if the first geometric matching condition (matching of 
the surface radius) ||2Jl is satisfied initially, then it is sufficient to demand that its coordinate time derivative 
^ is satisfied at all times. Secondly, with our choice of variables the second geometric matching condition 

is nothing more than an identity — with H22(l it has already been used to eliminate one variable. Thirdly, 
as will be shown in section^] equation H18|) . the time-component of the Lanczos equation, is in fact identically 
satisfied if all the other matching conditions are satisfied, the Einstein-field equations are valid on each side, 
and the surface-matter evolution is given by (|21|l . 

We conclude that the matching conditions are completely represented by ijSJ , the coordinate time deriva- 
tive of the first geometric matching condition, and the angular component of the Lanczos equation (|17(l 
together with an initial matching of the proper surface radius w = If. With our choice of variables these 
equations take the (surprisingly symmetric) form 



[L] = [uw + jwWu^ - 1] = (25) 
[wK\] = [uw' +jwVu'^ - 1] = -E, (26) 



where w = If is the proper surface radius of the spherical junction surface on each side and 

E ^ (27) 

quantifies the energy-content of the layer. To find a relation between K^g and L we square ()23(l and substitute 
from the square of (|24|l and obtain 

iwK\)^ ^L^ + a, (28) 

where a = w''^ — ilP . Versions of this equation have been given in [H] and For £' 7^ we can express 
Q in terms of L by using an algebraic identity as 

= [j-K'^fW-^'eV . (29) 
2\wK0h\ -2E ' ^ ' 



where 6 = a+ — a_. The explicit expression for L in terms of E takes the form 



Note that we find by differentiating H23() and H24|) with respect to u (taking w, w, and w' to be independent 
of u) the helpful relations 

jy/u^^^wK'g and ^^4^^^ - L, (31) 
ou ou 

which are valid on each side of the junction surface. 

The geometric matching condition (|24|l can be solved for u± and j± in terms of the time derivative of 

the surface radius L and the extrinsic curvature component g 

wL — w'wK^e 



—a 

We note that differentiating H32|) with respect to u and using H31|) yields 

w' L — wwK' 



(32) 



J 



vV^^^^iii^i^^^, (33) 



what also determines the sign of j and hence the radial direction of motion of the junction surface for each 
side. 
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2.5 The 'no surface-layer' case 

If all tangential components of the extrinsic curvature are continuous at the junction surface then it follows 
from l|26() that the surface-energy density ps vanishes. In this case the junction hypersurface is called a 
boundary-surface. It is an immediate consequence of (|28|l that 

h — — ^ — =y) ^ b = [a\ — [w — w J = 0, 

and 6 = ^ S = or = wK%/2 = ~wK%/2. 

Feasible solutions need to satisfy the two geometric matching conditions © and the derivative of 
the first matching condition (jHl and the matching of the extrinsic curvature l|17|l . Recognizing the similar 
structure of © and l(T7|) we form two new equivalent equations by adding and subtracting the two equations. 
The result reads 

[N{1 - al){w' - w)] = [N{1 + al){w' + w)] = 0. 

Using the factorized form of the second geometric matching condition Q and defining q = N{1 — al) this 
becomes 

[q{w' - w)] = [- {w' + w)] = 0. (34) 

If both w' — w and w' + w vanish separately on both sides, then the system becomes an identity and the 
junction surface motion remains undefined. Let us assume now that this is not the case. 

We note that for [w'] = [w] = the system is solved for any q+ = g_ . Hence, if the angular component 
of the metric and its first order proper time and radial derivatives are continuous, then the junction surface 
motion does not follow from the matching conditions. In particular, this is the case for the trivial matching 
of two identical space-times, were we have an 'imaginary junction surface', which could be placed anywhere. 

Let us from now on assume that at least one of the proper derivatives of w is not continuous at the 
junction surface. 

If w' ±w is zero on one side, it has to be zero on the other side too (otherwise no matching is possible) 
and one of the equations 1)34(1 is identically satisfied. 

However, if all w'j_ ± w± are non-zero then the condition 6 = also guarantees that the two linear 
equations H34|l are linearly dependent. The motion of the boundary surface is then described by the four 
equations (if w' ±w is zero then one of the options in the first equation is undefined, but the other option 
is then still valid) 

iVi (1 — aiZl) w' — w', + w+ „ , „ „ , db , , 

N^{l-a^_E) w'_^-w+ w'_+w- ±y ± ±' dt ' ^ ' 

2[Nw'{{w')' + w"al} - wN{w + {wYal}]. (36) 



where 



db 
di 



2.6 Expansion for small surface-energy densities 

As will be seen in the numerical examples given later, in many cases the surface-energy density (and hence 
E) approaches zero at some finite coordinate time. In the case of & 7^ the dynamic quantities can be 
approximated by a series expansion in terms of E. We start by re-writing the exact expression for the 
extrinsic curvature 1)29(1 as 

n b E 



It follows then from ((30(1 that 



2E 



where 0{E^) represents terms of the order E^ or smaller. Furthermore, from ((32(1 we find for u = N/k the 
expansion 

w'^b + sgii{L)w±\b\ 1 sgn{L)w±{a+ + a-)/\b\ T 3 
"± = 2^^ E + 2^1 ^ + )• 
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As the surface-energy density approaches zero the lapse functions (given by m+ and u-) diverge and the 
proper speed of the junction surface approaches the speed of light quadratically since 

\a±i±\ = 1 - ^Tir^n^-iu^ + ^ ^ • 

w^b + sgn(Ljw± \b\ 

To examine if and how E approaches zero we finally expand the evolution equation l|21ll 

Generally the first term diverges as we approach the speed of light — for example in the case of a perfect 
fluid one finds 

We conclude that if the energy-momentum contribution \T^^u^n'^] has a sign opposite to E, then E acceler- 
ates towards zero. In many cases E will reach zero at some finite coordinate time t^. Close to this point and 
assuming that the time dependence of all other terms is negligible we have £^ oc (<o — tY^^ . This implies that 
the lapse functions are integrable and the junction surface reaches the speed of light (on each side) within 
a finite proper time. Here our formalism breaks down and one would need a separate treatment of these 
singular cases. We want to speculate here that at these points the junction surface turns spacelike. 

On the other hand, if the sign of \T^^u^n^] is the same as the sign of E then E cannot get arbitrarily 
close to zero. In some cases (see figure^ E will oscillate around some value (which is itself time dependent). 
Even in these cases we can encounter divergencies resulting from diverging a± and b. This can lead to non- 
integrable lapse functions - from each side the junction seems to exist forever, but an observer who moves 
along the junction encounters a singular point after a finite time. At this point the surface energy density is 
zero and again the formalism breaks down. 

It should be noted that in the case of a perfect fluid on both sides of the junction the sign of the stress- 
energy contribution \T^i,u^-ny] depends on the energy density p and pressure p on both sides, i.e., on the 
equations of state. Hence whether a particular junction reaches the speed of light within a finite time or not 
might depend on the equation of state on each side. In section El we give a numerical example for such a 
case. 

Because points on the junction surface are not causally connected a spacelike junction surface has a 
very different physical interpretation . In such cases the junction surface cannot be treated as an 'evolving 
system' on its own, but rather as some kind of (spacelike) transition surface which is generated by the physics 
underlying the cosmological model. 

Usually a timelike junction surface is used to model the time evolution of a spatially localized inho- 
mogeneity. If a junction surface turns spacelike a breakdown in the thin wall approximation must have 
occurred. 

3 Necessary and sufficient conditions for a possible matching 
3.1 Demanding real solutions for L 

From (|30|l we find with > a necessary condition for the existence of solutions which restricts the allowed 
values for E, such that 

E^ -2E^{a++a^) + b^ >Q. (37) 

The roots of the quadratic polynomial (in E"^) on the left-hand side are given by a+ + a_ ± 2^a+a-. For 
a+a_ < or a+ < 0, a_ < all values of E are feasible. Hence it is only for a+ > 0, a_ > that restrictions 
on E arise, which take the form 

< |£^| < Vo^l or ^+^/^_<\E\. (38) 

The shape of the forbidden region takes a particular simple form in the _E — plane, which is illustrated in 
figure 13 

The two disjoint regions allowed for \E\ given by H38|l are easily distinguished by E^ < \b\ and E'^ > 
\b\. Furthermore, from ^ we find sgn{K'^g) = sgn(i4:^g) = -sgn{b/E) for E^ < \b\ and sgn(i4:%) = 
sgn(i?),sgn(if«g) = -sgn(i?) for E^ > \b\. 
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Figure 2: For the case a+, a_ > this figure shows the region in the — -^-plane with z =^ ■\/2(o+ + a_) 
where there is no solution for the junction surface motion. The dashed hues show the curves of constant 6, 
the circles are the lines of constant L/z, and the shaded region is not feasible since by construction 2|6| < . 



3.2 Proper time relations 

By setting fc = 1 it follows from Q the condition 

7V = u>+l, (39) 

i.e., on each side of the junction surface proper time must proceed faster (with respect to the time coordinate) 
than on the junction surface. To investigate the resulting constraints on the surface energy density we start 
by noting that 

\ 1 1 T^Q \ I I . -r I f > for a > . , „^ 

I^^ljio for aio ' (4^^) 

which can be easily verified by squaring and substituting from (I28II . Substituting (|32|l the inequality (|39(l 
takes the form 

(w'{wK\) - a) - wL ^ ^ 
a ~ 

It follows from H4n|l that for a > we need w'{wK^g) — a > 0, while for a < 

wL > 0. (41) 

Let us first consider the case a > 0. Using (|29|l the condition becomes w' {b + aE'^) /{—2E) — a > 0, where 
a = ±1 corresponding to the outside (+) and inside (— ) case. In the following let us use the convention that 
if a, w refer to the quantities on one side, then a,, refer to the quantities on the other side of the junction. 
By setting x = sgn(b)E/ \/\b\, e± = ±sgn(&) = sgn(a — a*), and s± = — we bring the inequality in 

the form 

- + ^-4??"f"^;n- (42) 

X [> 2s± for wit < 

The allowed ranges for x are illustrated in figure |2| 

The case of e > Let us note that if e+ > then e_ < and vice versa (e* < 0). The sign of the surface 
energy density is now determined by 

sgn{ps) = -sgn{b)sgn{w'). (43) 

Furthermore, if \s\ = \a/{w' ^/\b\)\ < 1 then no restrictions are placed on (but further restrictions could 
come from a* > 0). For \s\ > 1 the allowed range can be found by setting = and hence \l/x + x\ = 
2cosh(z). Using cosh~^ \s\ — ln(|s| + Vs^ — 1) one obtains the ranges 



< |a;| < \s\ - _ I |s| + - 1 < 
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(a) Allowed ranges for i? if e > 0. 



(b) Allowed ranges for i? if e < 0. 



Figure 3: Visualization of the inequality H42|l for e > (figure 3(a) I and e < (figure 3(b) I. The lower and 
upper horizontal line correspond to w'j^ > and w'j. < 0, respectively. Here s± = —a±/{w'j.^/\b\) are the 
terms on the right-hand side of (|42ll . The allowed ranges are indicated by a bold line. The end-points are 
given by {A and B are only defined for 
\s\ + Vs^ + 1. 



>1)A= \s\ - V^^~T, B = \s\ + V^^~l, C = -\s\ + V^^TT, D = 



The case of e < Similarly to the last case, if e < then e* > 0. Hence each case will occur once at the 
junction. A similar procedure as above yields the restrictions 

|s| — -f 1 < sgn(u'')a; < or \s\ + \J s'^ + 1 < sgn{w')x. 

The case of a > on both sides If a is positive on both sides then x can only take values which lay in 
the intersection of the allowed ranges on each side. The allowed intervals differ, depending on the signs of 
w'j^ and 'w'_. All possible cases are shown in tabled We note that in every case the allowed values for E 
have the same sign. 

As a particularly important case (for the matching of FLRW models) and an illustrative example we 
evaluate the restrictions on the surface energy density for a+, a_ > and w^, w'_ > 0. In this case x has to 
be negative and hence sgn(ps) = sgn(_E) = — sgn(6). If one assumes on physical grounds that ps should be 
positive then no matching will be possible if b is positive. 

The case of a < If a < on one side of the junction surface then (|41|l implies with sgn(L) = sgn(w) 
the sign for L, the coordinate time derivative of the surface radius, which was left undefined in H30|l . For 
the case that a+ and a_ are negative this condition must hold on both sides and hence a matching is only 
possible if 

sgn(w+) — sgn(ii;_) for a+,a_ < 0. (44) 

4 The time-component of the Lanczos equation 

So far we have only considered matching of the metric and of the angular components of the extrinsic 
curvature of the junction surface. The remaining matching condition comes from the time-component of the 
extrinsic curvature H18|l. which contains a second order time derivative of the junction coordinate radius, or 
equivalently a first-order time derivative of the lapse function. 
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b>0 

=^ e+ = = -1 


6 < 
e+ = -1, e- = +1 


w'j^ > 
w'_ > 


sgn(ps) = -1 
y|6|max(-A+,-C_) < E < for |s+| > 1 
— \/ 6 C_ < i? < for |s+ < 1 


sgn(p3) = +1 
< E < ^y\b\mm{A^,C+) for \s-\ > 1 
< -E < \/|b|C+ for |s- < 1 


> 
w'_ <0 


sgn(ps) = -1 
E< y|ft|min(-B+,-D_) for > 1 
E < — -^/|b|_D_ for s+l < 1 


sgn(ps) = -1 
E < i/|b[mm(-_B_,~_D+) for > 1 
E < —^y\b\D+ for s_ < 1 


w'j^ < 
w'_ > 


sgn(ps) = +1 
y/]b\max(B+,D-) < E for > 1 

/ 1 L 1 n ^ 77^ f 1 ^ 1 ^ 1 

\/\o\JJ— < £/ tor s+l < 1 


sgn(p3) = +1 
y|6|max(B_,D+) < E for > 1 
\/|o|-L'-j- < -C/ tor \S— \ < t 


w'_^ <Q 
w'_ <0 


sgn(ps) = +1 
<E < \/|fe|min{A+,C-) for > 1 
< E < \/W\C~ for < 1 


sgn(p3) = -1 
y/|6|max{-yl_,-C+) < £ < for |s_| > 1 
-v/WC+ < -B < for |s_ 1 < 1 



Table 1: For a+ > and a_ > this table shows the aUowed region for E for all possible combinations of 
sgn(?i;^) andsgn(w^). Here ^± = |s±|-y^4 - 1 < 1;B± = |s±| + y^4 - 1 > 1;C± = -\s±\ + ^sl + 1 < 

l-D± = \s±\ + .^4 + 1 > 1. and s± = -a±/{w'±^\). 



Rewriting the time-component of the extrinsic curvature in terms of our variable u — N/ k yields 



(45) 



where the factor N' /N is independent of the junction surface motion. Taking the coordinate time derivative 
of if^ and using (^5)) and ll^ we obtain 



j±LdU± 1 db dE 1 g 

^=^^ 2Edi He""^^' ^' 



4-1 



(46) 



where (note that {d/dt)f{t,a{t)R) = uf + jy/u'^ - If) 



,f dw' 



-dw 



z = u V j\Ju^ - 1 = ju\Ju^ - \{w" + ?«} + u^(w')' + iu^ - \){w)' . 

dt dt 

Differentiating lf77|) and using lf?T)) with = 2L/w yields 

dE 



dt 



which expresses the coordinate time derivative of E. Substituting for the first term in (|45|l allows us to 
evaluate the remaining junction condition [K^'^u^u^,] — —k{ps/2 +Ps)- The terms containing the surface 
pressure and density cancel each other and we obtain 



= -wlT^'^u^n,] + [z]~L 



i\/ V? ~ 1 



L 



(47) 



The first term can be expressed in terms of the Einstein-tensor with respect to the original metric as 



where the relevant components of the Einstein-tensor are given by 



Gtt 

Gtr 



W 

2/2 



,, / 1 — a 
w + w- H — - — 
/ 2w 



,N' I- a 
w ~ w — — H — - — 
w \ N 2w 

2NI ( N' 
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and L is given by H24|) . Substituting into 147|) and using the relation 

/ N' 
{w) -{w) jw - —W 

shows that (|47|l is an identity, satisfied for all sphericaUy symmetric junctions between solutions of the 
Einstein-field equations if the geometric matching conditions O and l@J) together with the angular com- 
ponent of the Lanczos equation (|26|l are satisfied. While it was well-known that for certain cases the 
time-component of the Lanczos equation is identically satisfied (e.g. |H1 it seems to be a new result for 
the generic spherically symmetric case. 

It was suggested that for the matching of FLRW models the time-component of the Lanczos equation 
determines the pressure 0. In light of the above result this cannot be the case and one needs to supplement 
the model with an equation of state for the surface-matter. 



5 Matching of FLRW sections 

We want to turn our attention now to the special case of the matching of two distinct FLRW regions. Such 
junctions are encountered in cosmological models which approximate universes containing many FLRW 
domains (multidomain universes). The most prominent example is Linde's Chaotic Inflation scenario ]5121- 
Junctions of this type have been studied in ^ JJ^ . Our treatment will serve as an illustration for the 
introduced method and as a source for numerical examples. Here it is not our aim to investigate physical 
processes which could lead the to creation of a "bubble" and we refer the interested reader to the vast 
literature (see, e.g., [Hl^^ElElEI)- Instead we want to focus on the generic geometrical and mathematical 
aspects. 



5.1 FLRW models and their parametrization 

The metric of FLRW models can be written in the form 

ds2 = -N^{t)dt^ + l^{t){dr^ + f{r)dn^}, (48) 

where l{t) is the scale factor, N{t) the so-called lapse function, dfl the line-element on the two-dimensional 
unit-sphere, and 

{sin(r) for closed models 
r for fiat models 

sinh(r) for open models 

Note that the FLRW metric H48|) has the same form as the general metric for spherical symmetric spaces 
(UJ, but with I' = 0, N'/N = and / = 0. 

The evolution of FLRW models is described by the Friedmann equation - the dynamic part of the 
Einstein-Field equations - 

/ .\ 2 

np + A C 



3 /2' 



(49) 



where C — 0, +1,— 1 for fiat, closed, and open models, respectively, and a dot indicates the derivative with 
respect to proper time t, i.e., I' = jf^- The matter is described by an energy- momentum tensor of perfect 
fiuid type. The unit tangent vectors to the fiuid flow lines are given by^ 

and the energy-momentum tensor describing the comoving perfect fluid takes the form T^'^ — {p + p)v^v'^ + 
pg^^ , where p is the energy density and p the pressure. The matter evolution is then described by the 
energy-conservation equation p+ [p -\- p)3H = 0, where H = l/l is the Hubble parameter. 

We restrict ourself to models with a 7-law equation of state, i.e., models in which energy density p and 
pressure p are related by"^ 

p=(7-l)p 7 e (2/3,2]. 

^In the new coordinate system which is continuous at the junction the components take the form i)f = -^<5^ — 
^The case 7 = gives an effective cosmological constant. We can exclude this case here because the cosmological constant 
is included separately. 
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In this case x-y ^pl^^ is a constant of motion. This aUows us to ehminate the energy density p from the 
Fricdmann equation (|49|l . so that the evolution of the scale factor I is described in terms of the constants of 
motion by 

H^^X-^l-^' + l^-ji- (50) 

Junctions are often used as models for transition regions between two almost-FLRW regions. Underlying 
is the assumption that an initially small transition region remains small. 

5.2 Comoving junction surface 

Let us first investigate whether there could be a comoving junction surface for two FLRW models with 
7-equation of state. This case is uniquely identified by a± = 0. 

From the geometric matching condition (specialized to the FLRW metric) Q we find 

/+("+) 

where the first factor is now time-independent. 

Let us assume we are given a solution to the Friedmann equation with J — ^-,x ~ X-; ^ ~ ^-i ^^'^ 
C = The question is now, whether there are constants 7+, A+, X'y^ and ^+ such that ^+ = A/_ is a solution 

for constant A > 0. Both solutions would have the same Hubble parameter iJ =' ^- = Substituting in 
each case from the Friedmann equation H5U|I gives 



-37- ,1a _ ^ 7-37+ -37+ , 1 A C+ 1 



which has to be satisfied for all values of L. On both sides all terms contain different powers of / (note that 
—37 € (—2,6]). In order that both sides contain the same powers of / we need 7_ = 7+. Comparing the 
coefficients gives then A = 1,C+ = C-i — ^-i ^^^d X7+ — Xi-- Hence the solutions are identical. 

We conclude that if the inside and outside of the bubble are evolving according to the Friedmann equation 
(gni with a 7-law equation of state then no non-trivial comoving junction is possible. We note that this 
result follows alone from the geometric matching condition (0} - it does not depend on the presence of a 
surface layer. 

5.3 Matching of FLRW regions with surface-layer 

The FLRW metric (|48|) implies w — l{t)f{r) and by taking proper-time and radial derivatives we derive 
w — Hw and ^' = The component of the energy-momentum tensor which is needed to evaluate the 
surface-matter evolution according to l|21|l is easily found to be (for completeness we include k, which is set 
to unity) 

'N'^al{p + py 



fc2 



We proceed now by expressing all quantities related to the metric and its derivatives {a^h,w' etc.) in 

terms of FLRW model quantities. With [Af /dr)^ = 1 — C,p we obtain the expressions a 1 — ^^^"^ w^. 

First we want to examine which kind of bubbles could exist if there can only be a positive surface-energy 
density on the junction surface, i.e., E > Q. For reasonably small bubbles (such that the circimference 
increases with the radial coordinate) we have a+ > 0,a_ > 0,andw' > 0. We find from table^that in this 
case 

sgn(-B) = -sgn(6) = sgn([Kp + A]), 

and hence junctions are only possible if the inside FLRW region has a smaller total^ energy density than the 
outside region. This is illustrated in figure 

*The cosmological constant represents the vacuum energy density. 
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sign(E)=+l 1 


sign(E)=-l 


[Kp] 

1 



Figure 4: Points in the A — Kp-plane for which a matching requires a positive/negative surface-energy density. 
Only certain matchings are possible if there can only be a positive surface energy density. 



5.4 The closed to (inflating) open junction 

To understand the behaviour of the junction surface it is instructive to consider a particularly simple example 
for which the evolution equations are known. One such example is the junction between a non-inflating 
closed geometry with radiation inside (A_ = 0, C- = +1) and an inflating empty open geometry outside 

(x+ = o,c+ = -i). 

For these cases the Friedmann equation H5U|) is easily integrated and one finds the well-known solutions 




where r+ and r_ are the proper times along fluid flow-lines outside and inside, respectively. For the inside 
model the scale factor L grows until it reaches a maximum at t_ = x_, and then declines until it reaches 
zero at r_ = 2%. On the contrary the outside model expands exponentially forever. 

1 f / 




Figure 5: Evolution of the scale factor in a closed non-inflating FLRW model (L) and in an empty open 
inflating model 



Since for each time t the proper time measured along the junction surface is always less or equal to the 
proper time in the inside and outside model {u± > 1) it is clear that such a boundary can only exist for a 
finite proper time measured along the junction — the (timelike) junction surface must be 'terminated' at 
some time. 
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There are four possible solutions. Firstly, it is possible that the junction surface exists forever (in terms 
of the proper time) in the outer region while the proper time along the junction surface is bounded. In 
our formalism this corresponds to a non-integrable divergence in the lapse function 7V_|_ = for the outer 
region. 

Secondly, the junction surface can contract to a point such that the inner region is eliminated. This case 
is characterized by Q!+, a_ and w approaching zero at some finite time. 

Thirdly, the closed surface might detach from the open geometry — the birth of a child universe. In this 
case the radial coordinate for the closed geometry approaches tt, while and w vanish (see figure IBJ. 



Figure 6: The closed inside geometry can detach from the open outside geometry. This happens when 
Q!_ — + TT within a finite time. 

As a last option the outside region might be eliminated. In the case of a closed outside geometry this is 
surely a possible solution, but in the cases of flat and open outside geometries this generally requires that 
the junction surface turns spacelike^. 

In our formalism such a behaviour would yield a diverging, but integrable, lapse function. At the 
singularity we reach the speed of light and our formalism breaks down. Nevertheless, of the three options, to 
become super luminal, to continue at the speed of light, or to decelerate, the first one seems most convincing, 
also with view on the results of section 1^1)1 

Generally, one of these cases has to occur before we reach the singularity in the inside region as can be 
seen from the following argument. Let us assume that for physical reasons only positive surface energies are 
allowed. Since the outside geometry is open we have w'j^ > 0. As the closed inside geometry approaches the 
'big crunch' singularity the energy density grows without bound. Hence b = w^{kp- — A+) has to become 
positive at some stage during the contraction phase. However, from table ^ one can see that there is no 
solution possible with > and 5 > if a+ and a_ are positive. Let us note that when the inside region 
is contracting we have W- < and w+ > 0. Hence according to 1)44(1 a+ and a_ cannot be both negative. 
If a-|_ is positive (and a_ negative) then l|43|l implies that the surface energy density is negative, which is in 
contradiction with our assumption. If on the other hand a+ < then this implies L > and hence the proper 
surface radius would increase. This just helps driving a_ = 1 — Kp-W^ /3 closer to zero, which eventually has 
to turn negative due to the diverging energy density. Again we reach a point where no solution is possible 
without negative surface energies. 

Note that if one allows negative surface energies then the above argument shows that if the junction starts 
with a positive surface-energy density then at some point the junction must have a vanishing surface-energy 
density. 

Figure IHl shows the results of a numerical integration of this particular model. It appears as if the speed 
of light is reached within a finite time (integrable lapse functions) on both sides. This strongly suggests to 
us that the junction turned spacelike. Note that this happens even far before the inner closed region enters 
the contracting phase. 

With this example we want to emphasize that there are junctions which are possible initially, but which 
evolve to some singular point. Numerical studies have shown that this is rather common. 

5.5 Vacuum bubbles 

As a simple case vacuum bubbles (p-f p = on both sides) have frequently appeared in the literature ll.2.j . 
For these cases the surface energy-momentum tensor takes the form 5"^,^ = —psS^i, and it follows from ((21|l 
that ps must be constant 

^One might speculate that for a closed inside geometry there has to be a finite volume and hence the junction surface has 
to turn backwards in time and become timeUke again. 
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Assuming '12' Kp + A =const. on each side one can integrate the evolution equation for the angular 
metric component (|30|) to obtain 

w{t) = Wo cosh.{t / Wo) , (53) 
where t is the proper time along the junction and 

v/(4(k(/9+ + P-) + A++ A_) + 3kV2)2 „ 64(kp+ + A+){kp^ + A_) ' 

which agrees with the findings in |12l . 

However, this does not yet establish the actual motion of the junction surface since the relation of the 
proper time (along the junction) to the coordinate time is unknown — we need the lapse function which is 
given by H32fl . The proper time along fluid flow lines (which is proportional to coordinate time) is then given 
by / Ndt. 



6 Numerical Results 

A computer program has been written to integrate the evolution equations for several FLRW junctions 
numerically. To achieve better accuracy around the singularities a variable step-width was used. All exam- 
ples given here are for positive surface-energy densities. Cases with negative surface energy can easily be 
constructed by exchanging the inside and outside region. The graphs on the following pages illustrate the 
results and will be discussed one-by-one below. 

Open inside, inflating closed geometry outside Figure[3shows such an example. After some time the surface 
radius starts to diverge (note the logarithmic scaling) while E = Kwps/2 does not approach zero (hence close 
to the divergence ps oc 1/w). The proper times on both sides seem to be diverging, which is in agreement 
with the results from subsection 12.61 Note that the inner and outer regions have a rather unusual equation 
of state with 7+ = 0.7 and 7_ — 1.9 — in this case it is really this choice of the equations of state which 
makes the energy-momentum tensor contribution [T^^u^n'^] positive for small values of E (see subsection 

Figure IHl shows the evolution of the same initial situation, but with different equations of state (dust on 
both sides). This seems to change the sign of the energy-momentum contribution [T^^u'^n'^] for small values 
of E, which now approaches zero within a finite time. In fact, it can be verified that close to the singular 
point to we have as expected E cx y^to — t. As predicted in subsection 12.61 the lapse functions appear to 
be integrable and the proper times do not diverge. The junction surface seems to reach the speed of light 
within a finite time. 

Clearly, our formalism breaks down at this point. However, one could argue that after reaching the speed 
of light within a finite time, one should expect the junction to turn spacelike. 

Closed inside, inflating open outside Figure El shows such a situation with a radiation equation of state for 
the inside region (7. = 4/3). This is the example discussed above. 

7 Conclusion 

We developed a formalism for the treatment of timelike junctions between spherically symmetric solutions 
of the Einstein-field equation, which is based on the Lanczos equation and the Israel junction conditions. 
We introduce new coordinates such that two conditions are satisfied: Firstly, all coordinates are continuous 
at the junction surface, and secondly, the junction surface becomes a surface of constant 'radial' coordinate. 
In this approach the actual movement of the junction surface is absorbed into the metric, of which the 
transverse components are discontinuous at the junction surface. 

We evaluate the junction conditions and re-discover with (|28|l and H29|l well-known relations between the 
extrinsic curvatures, the surface layer energy density, and the rate of change of the surface radius of the 
junction surface. It should be pointed out that these results follow without using the time-component of the 
Lanczos equation. As it was shown in section 0] for all spherically symmetric cases this remaining equation 
is in fact an identity. This was known for special cases, but it appears to be a new result in this general 
form. 
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energy density ps hj E = Kuipsfi. 



Figure 7: Evolution of a junction between an open (inside) and an inflating closed (outside) geometry. (Parameters :7s = 1,^+ = +1,C- = — 1)X+ = 5,X- = 
2,7+ = 0.7, 7_ = 1.9, A+ = 2,A_ = 0; Initial values: = 0.031, a+ = 0.20272, a_ = 0.2,/+ = !,/_ = l,p+ = 0.597, p. = 0.239. ) Note that here u and w 
are plotted logarithmically. 
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(a) Coordinate radius of the junction surfa<;e (b) Proper speed of the junction surface on (c) Lapse function N = u on each side 

on each side each side. 




(d) Proper times on both sides (e) Surface radius {w). (f) E, which is related to the surface-matter 

energy density ps hy E = Kwps/2. 

Figure 8: Evolution of a junction between an open (inside) and an inflating closed (outside) geometry. The lapse functions are integrable and the speed 
of light is reached within a finite time. (Parameters :7s = = +1,C- = ~liX+ = 5,X- = 2,7+ = 1,7- = 1,A+ = 2,A_ = 0; Initial values: 

Ps = 0.031, a+ = 0.20272, a_ = 0.2, Z+ = l,l_ = l,p+ = 0.597, p_ = 0.239. ) Note that here u and w are plotted logarithmically. 




Figure 9: Evolution of a junction between a closed (inside) and an inflating open (outside) geometry. The inside has a radiation equation of state, while 
the outside is an inflating dust model. (Parameters: 7,, = 1,C+ = ~1)C- = +1)X+ = 0,X- = li7+ ~ li7- = 4/3,A_|_ = 5,A_ = 0; Initial values: 
Ps = 0.005, a+ = 0.19739, a_ = 0.2, Z+ = 1, Z_ = 1, = 0, p- = 0.119. ) Note that here u and w arc plotted logarithmically. 



The behaviour for small values of = Kwps/2 has been investigated. It was shown that for certain 
cases E is driven to zero within a finite coordinate and proper time. At such a point our formalism breaks 
down. Nevertheless, we want to speculate here that in such cases the junction really turns spacelike. This 
can be seen as an inadequacy of the thin wall formulation in such situations — a causal propagation of a 
discontinuity should not exceed the speed of light. We suggest that in such cases the spatial extent of the 
transition region is not negligible. 

The developed formalism gives us two sources for constraints on possible junctions. Firstly the time 
derivative of the surface radius is given by the quadratic equation H30() . Demanding that real solutions to 
this equation must exist directly restricts the possible values of the surface energy density for a particular 
junction (see figure EJ. Secondly, in our approach physical solutions must have a lapse function which is 
greater than or equal to unity. The resulting restrictions depend on the metric components on each side of 
the junction — they either determine the sign of the derivative of the proper surface radius, or they restrict 
the possible surface-energy densities (see figure ISJ. For the latter case the allowed ranges for E — nwps/'2, 
have been given explicitly. 

For the special case of junctions between FLRW models with 7-equation of state it was shown that alone 
on geometrical grounds there can be no comoving junction surface — whether with or without surface layer. 

As a particularly simple and well-known case, 'vacuum bubbles' were discussed, and the results agree with 
the literature. A particularly interesting model, the junction between an empty, open, inflating FLRW region 
outside and a radiation dominated closed FLRW model inside, has been investigated in more detail. The 
inside region re-collapses after some finite proper time and hence the junction surface has to be terminated. 
Besides a disappearance or a detachment of the closed inner region we suggest that the junction can turn 
spacelike — an effective disappearance of the outer region. 

This and other examples have been integrated numerically. It was observed that many models seem to 
reach the speed of light within a finite proper time, in accordance with the predictions from section 12.61 
Since a spacelike junction violates causality, we suggest that a breakdown in the thin-wall approximation 
must have occurred. 

Our results show that the thin- wall treatment of timelike junctions (without the presence of scalar fields) 
is on a mathematical sound level. Nevertheless, in many cases the junction surface reaches a singular point 
within a finite proper time. We believe that in these cases the thin-wall is not a physically acceptable 
approximation. 
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